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THE  CALCULATION  OF  INCLUSIVE  GAMMA  SPECTRA  AT 
HIGH  ENERGIES 

R.  Blutner,  S.  Frank,  K.  Hansgen,  Htim  Than,  H.  KSnig, 

J.  Kripfganz,  H.  Kuhnicke,  H.  Lemke,  H.  J.  MShring,  V.  Muller, 
G.  Ranft,  J.  Ranft  and  F.'  Thiel 

1 . Introduction 


Gamma  quanta  vdiich  are  produced  in  highly  energetic  pp 
or  v collisions  arise  almost  exclusively  from  the  decomposition 
of  the  -71^  mesons  which  are  produced  first: 


I — > y+y 

np  71^  "I"  **' 


(1) 


This  fact  makes  it  possible  to  find  the  differential  effective 
cross  section  of  the  Y quanta  produced  by  purely  kinematic 
calculations  from  a known  7i°  spectrum.  The  spectrum 

is  calculated  according  to  the  thermodynamic  model  [2,  4,  5, 

6,  7,  9,  10].  This  takes  both  of  the  essential  possibilities 
for  production  of  mesons  into  consideration; 

The  arises  through  direct  decomposition  of  a fireball. 

The  fireball  is  considered  as  a resonance.  This  resonance 
decomposes  into  7i°  and  j>- 

In  Section  2 we  calculate,  for  the  special  case  of  isotropic 

(1)  This  work  was  produced  voider  the  care  of  Professors  G.  Ranft 
and  J.  Ranft  by  a student  group  consisting  of  7 second  and 
third  year  students  (as  part  of  the  WPS),  2 graduates,  one 
research  student  and  one  aspirant.  At  the  Exhibition  of 
The  Karl  Marx  University,  Leipzig,  1972,  it  was  distinguished 
with  the  Prize  of  the  Karl  Marx  University,  1st  Stage. 

1 


PTD-ID(PS)T-n9nO-77 


71°  spectra  an  expression  for  the  differential  effective 
cross  section  of  the  Y quanta. 

For  the  general  case  of  anisotropic  distributions 
we  give  in  Section  3 a Monte  Carlo  method  with  which  the  Y 
spectra  can  be  calculated.  In  order  to  determine  the  Y 
spectra  with  good  statistics,  however,  the  Monte  Carlo 
program  requires  a substantial  computer  time.  Therefore, 
in  Sections  4 (direct  production  of  the  n'  ) and  5 
(production  of  the  71°  through  resonance)  the  spectrum  is 
in  part  analytically  calculated.  Our  results  are  reported 
in  Section  6 and  discussed  in  comparison  with  the  existing 
experimental  data. 

2.  The  Calculation  of  Y-Spectra  from  Isotropic  ■ Spectra 

Let  an  xms table  particle  of  mass  m be  generated  in 
a collision  between  particles  a and  b.  Let  the  momentum 
spectrum  of  this  particle  be  isotropic  in  the  center  of 
gravity  system;  that  is, 

■7^,  =/(l?l)  ■=/(/’)  ~ >n*^)  ■=/(«).  (2.1) 

Let  the  particle  of  mass  m decompose  into  particles 
1*  and  j with  masses  mj^  and  mj . We  seek  the  momentum 
spectrxm  1/,  (e)  of  the  particle  1 in  the  center  of  gravity 
system. 

Under  the  assumption  that  the  unstable  particle 
disintegrates  isotropically,  fM  can  be  determined  from 
purely  kinematic  considerations.  These  calculations  are 
presented  in  [3].  If,  instead  of  the  momentum  distributions 
fM  and  f(e)  one  introduces  the  energy  distributions 


W(£*)  = 

(u,(£)  = 47rpcf(e)  (2.2) 

[Throughout  the  document  1 = Italic  Z.] 
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then  the  result  is 


W,(£) 


m 

'^Pt 


* p** 

n J 


oi{e'*) 


(2.3) 


with 


• 

«*  = — r ± ^/’i) 

mf 

w**  + (m?  - /nh 

Cj  — £ 

Im* 

r-i  - r*  _ .f,2  _ t'"**  “ ('”i  + 

Pi  ^ — 

4m*^ 

These  relations,  however,  apply  only  if  the  rest  masses  of 
particles  1 and  j are  finite  (m,  + o,  jmj  + O). 

On  disintegration  of  the  tmstable  particle  of  mass  m 
into  particles  1 and  j with  vanishing  rest  mass  (e.  g., 

Y quanta,  a limit  transition  is  necessary,  m, = m ->  0 
This  limit  transition  gives,  for  the  limits,  instead  of  the 
expression  (2.4) 


(2.4) 


- - mjY] 


£_  = C + 


At 


(2.5) 


For  finite  collision  energy,  we  must  use  instead  of  e? co 
the  maximum  possible  energy  of  the  unstable  particle; 


£ 


(2.6) 
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The  Monte  Carlo  Procedure 


If  the  momentum  distribution  of  the  unstable  particle 
of  mass  m is  anisotropic,  then  (2.3)  must  be  replaced  by 
a considerably  more  complex  multiple  integral  which  can  be 
evalmted  only  by  means  of  drastic  approximations  (e.  g.,  [ll]) 
or  tedious  nvmierical  calculations.  In  this  case  it  is  simpler 
to  calculate  the  Y spectrum  using  a Monte  Carlo  procedure. 

In  the  Monte  Carlo  procedure  the  scattering  or  disin- 
tegration process  is  simulated  in  a computer.  Here  we  consider 
the  decomposition  =1° -» 2>'  . The  -'t°  spectra  can  be  specified 

in  arbitrary  form  (e.  g.,  from  an  empirical  formula,  or 
according  to  the  thermodynamic  model  [2,  4,  5,  6,  7,  9,  10]. 
Corresponding  to  the  selected  statistics,  from  this  one 
determines  the  number  of  mesons  which  vanish  in  a single 
phase  space  element.  In  the  rest  system  of  the  the 
disintegration  ji® 2y  is  isotropic  ( has  zero  spin).  The 
disintegration  angle  is  a free  parameter  and  in  the  ^tonte 
Carlo  procedure  it  is  occupied  by  random  numbers . Then  the 
disintegration  is  unambiguously  defined.  The  momentum  of  the 
Y quanta  is  transformed  back  to  the  starting  system. 
Corresponding  to  the  phase  space  decomposition,  the  ^quanta 
are  assigned  to  their  phase  space  elements.  Thus  one  obtains 
the  Y distribution  as  a histogram. 

4.  Y Spectra  According  to  the  Thermodynamic  Model;  the 
Contribution  of  Directly  Generated  Mesons 

For  the  thermodynamic  calculation  of  the  contribution 
to  the  Y spectra  resulting  from  the  disintegration  of  directly 
generated  mesons  we  need  the  function  /°(f,  H This 
function  designates  the  m ntum  distribution  of  the  Y qiaanta 
arising  from  directly  generated  .i"  mesons.  Here,  the  reference 
system  for  the  resting  system  is  a fireball  of  temperatvire  T. 
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According  to  the  thermodynamic  model  [2,  4,  5,  7,  10]  the 
spectrum  of  the  directly  generated  .■»*’  mesons  in  the  resting 
system  of  the  fireball  is  isotropic: 

Z*V  1 

- 8^3  • ^.V>  _ 1 • (4.1) 


In  the  following,  we  neglect  the  one  in  the  denominator. 

If  the  51°  disintegrates  isotropically  into  particles  of  mass 
mj^  and  m^,  then  an  elementary  calculation  of  the  integral  (2.3) 
gives 


/<(0  - - 


m*z*z,  V 


(Cf.  [3])(4.2) 


Here  m is  the  mass  of  the  unstable  particle  (in  our  case, 
of  the  5t°  meson);  z-,  and  z are  the  multiplicities  of  the 
particles  with  masses  mj^  and  m , respectively. eiand 
Pj^  are  given  in  Formula  (2.4).  The  Y spectrum^in  the  rest 
system  of  a fireball  arises  through  the  limit  transition 

/n,  — ntj  = 0. 


As  is  shown  in  Section  2,  in  (4.2)  we  must  use  for 
«-  and  «+  the  expressions  (2.5)  and  (2.6),  respectively: 


/;  c.  n - ^ r>  [(l  + . ■■  ■<'>  - t .-..vn  (4  3) 


5.  Y- Spectra  According  to  the  Thermodynamic  Model;  The 

Contribution  of  Mesons  Generated  by  the  Decomposition 
of  Continuous  Resonances 


5.1  Formulation  of  the  problem  and  transformation  of  the  integral 

We  need  the  fmction  /"(c,  r)  , which  is  the  analog  of  the 
function  /•(«,  T")  . But  here  the  mesons  are  generated  by 

the  decomposition  of  a continuous  resonance.  Thus,  we  consider 
the  disintegration  scheme  in  (1).  In  order  to  obtain  the 
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momentijm  spectrum  of  the 


mesons  from  (4.2),  we  set  there 


Hlj  = w, 

Wj  = m, 

jii*  = mass  of  the  resonance 


(5.1) 


The  momentum  spectrum  /*  (e,  T)  of  the  Y quanta  in  the  rest 
system  of  a fireball  arises  from  this  through  a further 
integration  of  (2.3)  and  (2.4).  Then  we  perform  the  limit 
transition  m,  0 : 


fy  - f 1 + — ) < 


^ITl 


(5.2) 


f 4- 


— C + = )fVl  + 

4c 


M"*  = 


(we,  ± Vw*  — s/ej  — /mJ) 


e.  = 


+ /nj  — /mJ 
' 


(5.3) 

(5.4) 


Next,  we  must  calculate  the  integral  (5.2). 
the  following  designations; 


We  introduce 


(5.5) 


Because  e_  > m.  and  e_  = m.,  only  if  e = m./2  , then  for  «+"»./2 
the  integrand  in  (5.5)  is  a continuous  function  in  the  range 
of  integration.  The  substitution 


= cos  /iz 

— = cos  ha 
in. 


(5.6) 
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yields 

+ ~cosh{z  + 

A further  substitution 

/_ ; cos  h{z  — d)  = y - 
: cos  /?  (z  + o)  = y + - 


gives 


h - f’"*’  7- dy. 
^ Jy,o-)  Vi -(]/y)\  ^ y) 


The  function  = cos /i(z(co)  + fl)has  the  form 


y,(a))  = --  -^  + - 1 /-4-  - ' • 

m,  m,  V '«x  V 

At  the  lower  limit  this  expression  simplifies  to 

\ ^ V + u falh  e > /nJ2 

m,  4e 

c w_  „ „ 

/,  ■)  = u H i;  falls  e < mjl 

4e 

/f  ) = u + i;  falls  e > iriJ2 

)*v->  m,  4e 

= — - V + u falls  £ < mJ2 


at 


(5.7) 


(5.8) 


(5.9) 


(5.10) 


(5.11) 


(5.12) 
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We  always  have 

y,(f.)  S 1. 

If  l(a,b)  designates  the  integral  (5.9)  with  the  lower  limit 
a and  the  upper  limit  b,  then  we  miist  solve  this  integral 
for  CO  > b > a > \ , 


5.2  Analytical  Solution 


If  the  lower  limit  in  (5.9)  is  large,  then  it  is  a 

t 

good  approximation  to  expand  the  root  in  the  denominator 
as  follows: 

1 


i'-^r 


1 + 


2y^ 


With  this,  (5.9)  gives 


By  virtue  of  the  relations 

j: 


d;; 

y 


-EA- 


(>n*  \ 


where  f,(A)  is  the  exponential  integral  and 

J ^ ! 1 ( 7-1  ' T J ■7'’  V’ 

it  is  possible  to  calculate  the  integral  (5.14): 

{a.  b 


, w* 

, M 

- £ ('^  A 

- 1 

\ 

1 + i 

47y 

4>-V 

\ 

(5.13) 


(5.14) 


(5.15) 


(5.16) 


(5.17) 


8 


One  can  easily  estimate  that  for  a > 1.4  the  error  arising 
from  the  approximation  (5.13)  is  less  than  15%. 


In  order  that  the  integrand  not  vary  excessively,  it  is 
desirable  to  divide  the  term  additively: 


{a,  b)  = w {a,  b)  + //>■  {a,  b) 


(5.20) 


— = /f 

''■‘“•‘>’11^^’  (5.22) 

Vx^  + 1 J 

(z)  = -ife-'’  d/ 

Here  v^Jo  is  the  error  function.  The  main  part 

of  w(a,b)  appears  in  (5.20)  if  a and  b are  near  1. 


5 . 4 Summary 


In  the  calculation  of  the  Y spectrm  in  the  rest  system 
of  the  fireball  (5.2)  the  integral  (5.9) 

/_  0'_(e-). ;'-(«+))  and  /+(}'+(«-).>'+(«+)) 

must  be  solved.  Depending  on  the  positions  of  the  lower  limit, 
a,  and  the  upper  limit,  b,  we  can  distinguish  3 cases: 

Case  1:  1.4~uusa<b 

Formula  (5.17)  is  used  with  an  error  less  than  15%. 

Case  2.  l^sa<bsuu=1.4 
JIs  {a,  b)  (5,22) 

(5.20)  is  used.  ^ here  must  be  calculated  numerically 
using  a Gaussian  integration  formula. 

If  the  upper  integration  limit  here  is  far  above  the  lower, 
then  one  can  arbitrarily  reduce  the  upper  limit  without  making 
a large  error.  But  the  error  occurring  in  the  numerical 
integration  is  smaller. 

Let  us  calculate 


A)m*IT{b  - a)  < D (D  specified  « 4) 

Js{o,b)  - ^^■{a,b)  /r,(a,b)  (5.23) 

b)  ib  ~ a)  > D 


Set  b = r/m*  D a and  calculate 

(a,  b)  - If  {a,  b)  f //^.  {a,  b).  (5.24) 

The  error  remains  smaller  than  e-®  joo%. 
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Case  3.  ls:aiuu=1.4ib 


This  case  is  reduced  to  cases  1 and  2: 


J {a,  b)  = Ik  {a,  uu)  + {uu,  b) 


(5.25) 


The  function  (5.2)  /y(e,T)  for  generated  by  resonance 
disintegration  is  computed  numerically  using  the  formulas  above 
by  means  of  a FORTRAN  program.  Using  the  program  SPUR  19] 
a contribution  to  the  differential  effective  cross  section 
for  Y production  is  calculated  from  this . The  other  contri- 
bution arises  from  the  function  /°  (£,/')  (Formula  (4.3)). 

Both  contributions,  added,  give  the  differential  momentum 
spectrum  for  Y production. 


6.  Results  and  Discussion 

The  version  of  the  thermodynamic  model  based  on  the 
"strong  bootstrap  solution"  (see  [6])  is  used  for  the  thermo- 
dynamic calculation  of  the  Y spectra.  The  free  parameters  of 
the  model  are  adapted  to  19  to  70  GeV/c  data  from  pp  collisions. 
We  assume  that  the  observed  Y spectra  can  be  completely  ex- 
plained from  the  decomposition  of  n®  mesons.  Our  results 
indicate  that  this  assumption  is  apparently  justified.  The 
thermodynamic  parameters  for  the  reaction  p p * a®  \- 
are  needed  for  the  calculation  of  the  y spectra.  As  there  are 
neither  thermodynamic  calculations  nor  experimental  data  for 
this  reaction,  we  use  in  place  of  the  spectra  the  thermo- 
dynamic n*  or  n spectra.  In  Figure  1,  the  y spectra 
calculated  with  the  Monte  Carlo  procedure  (Section  3)  in  the 
form  <i^A7(5;)dC  in  the  center  of  gravity  system  are  compared  with 
the  experimental  data  [1]  for  the  reaction  p+p  -^y  + --  at 
23.1  GeV.  For  the  calculation  of  the  thermodjmamic  Y spectra 
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it  is  assvimed  that  the  .n'  spectra  are  identical  with  the  51+  spectra. 
Figure  2 shows  the  comparison  of  the  same  data  with  the  direct 
thermodynamic  calculation  which  was  described  in  Sections  4 
and  5.  Here  it  was  assimied  that  the  spectra  are  identical 
with  the  spectra. 


Figure  1.  Comparison  of  the 
experimental  Y spectra  [ 1] 
in  the  form  Niw>dQ) cm 
p<„b-=  2^,1  Gev/c  with  the 
thermodynamically  calculated 
Y spectra.  The  Monte  Carlo 
method  was  used  for  the 
calculation.  We  used 
spectra  instead  of  the 
spectra. 


Figure  2.  Comparison  of  the 
experimental  Y spectra  [l] 
in  the  form  {d^Nicpen)cm  at 

= 23,1  GeV/c  with  the 

thermodynamically  calculated 
Y Spectra.  The  analytical 
method  was  used  for  the 
calculation.  We  used 
spectra  instead  of  the 
spectra. 
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Inclusive  Y spectra  p ^ p ♦>'  + •••  were  measured  in  an  ISR 
experiment  [8].  The  collision  energy  corresponded  to 
laboratory  momenta  of  up  to  1,500  GeV/c.  Figures  3 and  4 show 
the  comparison  of  these  data  in  the  form  d'^NldpoQ  in  the  center 
of  gravity  system  with  the  analytically  calculated  Y spectra. 
It  is  likewise  assumed  that  the  spectra  are  identical  with 
the  71-  spectra.  In  Figure  5,  one  of  the  measured  curves 
at  E =52.7  GeV  is  compared  with  the  result  of  the  Monte 
Carlo  method;  as  in  Figure  1 the  spectra  are  assumed  to 
be  identical  with  the  spectra.  In  Figure  6,  the  data  at 
= 52.7  GeV  in  the  foirm  i^ii^Njd^p)  at  /^n  = o are  compared 
with  the  results  of  the  Monte  Carlo  calculation. 

We  could  establish  that  both  methods  for  calculating 
the  inclusive  Y spectra  according  to  the  thermodynamic  model 
yield  consistent  results,  as  can  be  seen  from  the  calculated 
spectra  at  23.1  GeV  in  Figures  1 and  2,  and  at  1,500  GeV 
in  Figures  4 and  5.  We  find  that  the  experimental  data  at 
23.1  GeV  approximately  coincide  with  the  calculated  values. 

The  deviations  may  be  explained  from  the  fact  that  we  used 
the  thermodynamic  ti'  or  spectra  instead  of  the 
spectra.  At  large  angles,  the  major  contribution  to  the 
spectra  comes  from  the  directly  produced  (Section  4). 

At  smaller  angles  the  contribution  from  71”  produced  in 
continuous  resonances,  calculated  in  Section  5,  increases. 

The  thermodynamic  calculation  of  the  Y spectra  at  1,000 
and  1,500  GeV  and  the  comparison  with  the  ISR  data  is 
significant  for  the  following  reasons: 

Conclusions  about  the  ;i"  spectra  at  these  high  energies 
can  be  drawn  from  the  Y spectra. 

Theyspectra  are  measured  up  to  angles  of  90®  in  the 
center  of  gravity  system.  The  model  is  compared  with 
experimental  data  for  the  first  time  at  energies  up  to 
1,500  GeV  in  the  central  region. 
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Figure  3.  Comparison  of  the  Y spectra  measured  by  ISR  [8]^ 

in  the  form  {d'‘Ni8peP:  cm  at  = 44.7  GeV  with 

the  thermodynamically  calculated  Y spectra. ^ The 
analytical  method  was  used  for  the  calculation. 

We  used  j*'  spectra  instead  of  the  spectra. 


Figure  4.  Comparison  of  the  inclusive  Y spectra  measured  by 
ISR  [8]  in  the  form  {d'‘Ni8iKO)  at  E = 52.7  GeV  for 
a scattering  angle  of  a™  = 10.2®, 90°  with  the 
thermod3mamically  calculated  Y spectra.  The 
analytical  method  was  used  for  the  calculation. 

We  used  t instead  of  the  i®  spectra. 
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Figure  5.  Comparison  of  the  inclusive  Y spectra  [8]  measured 
by  ISR  in  the  form  for  ® ’^.2'  at  an  energy 

= 52.7  GeV  with  the  thermodynamically  calculated 

Y spectra.  The  Monte  Carlo  method  was  used  for 
the  calculation.  We  used  spectra  instead  of  the 
spectra. 


Figure  6.  Comparison  of  the  inclusive  Y spectra  [8]  measured  by 
ISR  in  the  form  £:•(()*  A'/dVj  as  a function  of  pI  for/’.. 
at  an  energy  of  Ecm  = 52.7  GeV  with  the  thermodynamic- 
ally calculated  Y spectra.  The  Monte  Carlo  method  was 
used  for  the  calculation.  We  used  instead  of  .v 
spectra. 
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The  thermodynamic  and  spectra  are  compared  at 

energies  of  19  to  70  GeV  with  the  data.  It  is  shown  in  [6] 
that,  for  small  angles  and  energies  of  500  to  1,500  GeV,  the 
extrapolated  thermodynamic  spectra,  including  their  absolute 
standardization,  approximately  agree  with  the  ISR  data  for 
inclusive  or  production.  In  our  work,  likewise,  the 
Y spectra,  including  their  absolute  standardization,  are 
predicted  in  the  entire  kinematically  possible  range  and 
compared  with  experimental  data  in  a significantly  larger 
angular  range  than  in  the  case  of  the  -z*  spectra. 

Figures  3 to  6 show  approximate  agreement  with  the 
experiment.  It  follows  from  the  agreement  obtained  that 
for  these  energies  and  the  angles  considered  the  differences 
between  the  n*-,  n-  and  spectra  are  apparently  slight. 

The  thermod5mamic  spectra  show  "scaling"  behavior  and  a 
non-vanishing  central  plateau  [6].  It  follows  from  the 
comparison  with  the  Y data  that  the  spectra  in  the 
fragmentation  and  pionization  region  apparently  also  show 
at  least  approximate  "scaling"  behavior.  The  deviations 
which  occur  between  the  thermodynamic  Y spectra  and  the 
experiment  we  explain  by  the  following  facts ; 

1.  Use  of  ^"^or  ji"  spectra  instead  of  spectra. 

2.  Extrapolation  of  the  thermodynamic  spectra  from 

lower  energies  or  small  angles  to  energies  of  1,500  GeV 
and  the  entire  angular  range  also,  in  the  case  of  the 
and  spectra  [6]  leads  to  differences  of  the  same 
order  of  magnitude  from  the  experimental  ISR  data. 

Exact  matching  of  the  spectra  with  the  ISR  energies 
was  not  attempted  in  any  of  the  cases . 


We  thank  Mr.  Magnus  for  support  in  programming;  Mr. 
Ilgenfritz  and  Mr.  Hatthaus  for  useful  discussions.  Not 
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least,  we  thank  Mrs.  Eichel  for  the  carefully  performed 
copy  of  the  manuscript. 

Sxjmmary 

In  high-energy  collisions,  Y quanta  arise  from  the 
disintegration  of  primarily  produced  mesons.  The  inclusive 
momentum  distributions  of  the  Y spectra  are  computed  through 
a transformation  of  the  spectra  by  means  of  the  thermo- 
d5niamic  model.  The  work  was  stimulated  by  a group  of 
Soviet  physicists  at  the  VIK  bubna,  who  performed  the 
corresponding  experiment  using  a propane  bubble  chamber  at 
the  76  GeV  accelerator  at  Serpuchow.  The  results  of  this 
experiment  should  be  used  for  comparison. 
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